Heat kernel of integrable billiards in a magnetic field. 



We present analytical methods to calculate the magnetic response of non- 
interacting electrons constrained to a domain with boundaries and submitted 
to a uniform magnetic field. Two different methods of calculation are consid- 
ered - one involving the large energy asymptotic expansion of the resolvent 
(Stewartson-Waechter method) is applicable to the case of separable systems, 
and another based on the small time asymptotic behaviour of the heat ker- 
nel (Balian-Bloch method). Both methods are in agreement with each other 
but differ from the result obtained previously by Robnik. Finally, the Balian- 
Bloch multiple scattering expansion is studied and the extension of our results 
to other geometries is discussed. 
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I. INTRODUCTION. 



The aim of this work is to present analytical methods to calculate the magnetic response 
of non-interacting electrons constrained to a domain with boundaries and submitted to a 
uniform magnetic field. 

Historically, this problem traces back to the Bohr-van Leeuwen theorem stating the 
absence of classical orbital magnetism due to the exact cancelation between the bulk and 
edge magnetizations Later on, Landau using a quantum approach did show the existence 
of a finite magnetization [Q]. An extension to finite systems was provided by Teller who 
showed that the Landau magnetization results from an almost cancelation between the bulk 
and edge contributions. In this work, we shall concentrate on the problem of non-interacting 
electrons using the semi-infinite plane as a paradigm (our methods can be readily extended 
to other integrable systems - for example a disc). 

Whereas it is important to know the spectrum with a sufficient precision in order to 
describe low-temperature and high magnetic field phenomena (such as the Integer Quantum 
Hall Effect), the high temperature or weak magnetic field response like for instance the 
orbital diamagnetism is determined using smoothed spectral quantities. We shall describe 
them here by defining and calculating the heat kernel or equivalently its Laplace transform. 
The small time asymptotic expansion of the heat kernel is simply related to the smooth 
part of the density of states (the Weyl expansion) |^ and to smoothed thermodynamical 
quantities like the magnetization This asymptotic expansion of the heat kernel for the 
semi-infinite plane gives the perimeter correction to the Landau diamagnetism as noted by 
Robnik . We shall expand the heat kernel in two different ways obtaining the same results 
which differ however from those obtained by Robnik. We subsequently compare our results 
with those obtained by using the Balian-Bloch expansion. We show (in agreement with 
more recent works ||^) that the Balian-Bloch method is not very convenient in the presence of 
a magnetic field due to the fact that several terms in the multiple scattering expansion have 
to be included in order to obtain just one term in the asymptotic series of the heat kernel. 
Another asymptotic method we use, following Stewartson and Waechter 0, is powerful 
enough to give in principle all the terms in the asymptotic series for the semi-infinite plane 
in a magnetic field. We argue that all the coefficients in these series are universal in a sense 
that we shall precise. 



II. GENERALITIES ON THE RESOLVENT AND THE HEAT KERNEL. 

The resolvent G{E) of a bounded system is by definition the Laplace transform of the 
heat kernel P{t) = Tre~^^* = X^n ©"S'^"*, i.e. formally: 

where the sum runs over all eigenstates n of the Hamiltonian H of the system. This function 
was extensively studied for the problem of the Laplacian on manifolds with boundaries 
1^,11]. As is well-known, the small t behaviour of the heat kernel of a two-dimensional 



system is given by Weyl's formula: P{t) = ytt + 0{t ■2), where S is the area of the system; 
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therefore its Laplace transform is actually not well-defined. A natural way to regularize it 
is to subtract this leading small t behaviour to P{t) and then take the Laplace transform: 



oo 



g{E) = \j dt{P{t) - ^)e-^^*, Re(E) > 0. (2) 



An expression for g{E) valid also for Re(E) < was given by Berry and Howls [Q: 



g(s = Ve) = lim 



M-^ 1 mS r En 



n<N 



(3) 



It is easy to show that g{s) is related to the density and to the integrated density of states 
by: 

d{E) = ^ 5 - E J = ^ - i hm Im g(i v^E^) (4) 



n 
E 



C{E) 



where the contour C{E) in the complex plane encloses the segment [0, E\ of the real axis (so 
that it contains inside a finite number of poles of g{i^/e)). Using (|^) and the fact that P{t) 
is the Laplace transform of d{E), one can show that the two definitions (0) and (|^) coincide 
for Re(E) > 0. 

Given one of the three asymptotic expansions for small t or large s, E, of P{t), g{s) or 
d{E), one immediately gets the two others. For example, assuming the Weyl expansion of 
the resolvent is known: 



9is) ~ E ^> ^ ^ ^> |Arg(^)| < - - A, A > 0, (6) 



r=l 



we obtain: 

1 ^ (-!)'■ ^ _ 

m-^,^t^^il)i'\t^,. (8) 

The equivalence between (|]), (|^) and (||) is a consequence of equation (^), the generalized 
Watson's lemma, and of its reciprocal [0. 

As a simple example, the heat kernel and the resolvent associated to the Landau spectrum 
of a particle moving on a two-dimensional plane in a perpendicular field are easily found to 
be 

Pooit) = (9) 



2 sinh f 



and 
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9oo{s) = ^{H-u-\)-^{-u)), (10) 

if the degeneracy of the levels is N<^ = SB/^o and ^ = ~^~\ {i'i.z) = is the digamma 
function). Their asymptotic expansions are: 

doo{E) ~ — -5- (no more terms). (11) 

The last result is obvious because smoothing out the regular step-like integrated density of 
states over the energies a linear function with no higher power-law corrections is obtained. 
Because for any two-dimensional system of surface area S in magnetic field, P{t) has same 
leading small t behaviour as Poo{t), another way to regularize the resolvent is to subtract 
Poo{t) from P{t) and then take the Laplace transform. Such a regularization gives g{s) — 
Qoois), and can be calculated in the following way. 

Let us define the Wick-rotated retarded Green's functions G^it; r, r') and G^it; r, r') by 
G+{t; r, r') = G+it; r, r') = if t < 0, and 

{H + ri^)G^it- r, r') = h5{t)5{v - r'), (12) 

with the additional Dirichlet condition for G^{t\ r, r') on the boundary OS of the system: 

G+(t;r,r') = for r G dS. (13) 

The Laplace transforms of these Green's functions, G{E\ r, r') and Goo{E] r, r'), are defined 
by the same boundary condition ([T3|) and the equation: 



{H + E)G{E;r,r') = 5{r-r'). (14) 
Then using and the definition of the heat kernel, we have: 

g{s) - g^{s) = J dh{ G{s^; r, r) - G^s'; r, r) ). (15) 



III. EXPRESSION OF THE RESOLVENT FOR THE SEMI-INFINITE PLANE. 



We propose to apply formula ([T5|) to extend for the magnetic case the method developed 
by Stewartson and Waechter 0. We consider the following setup. A spinless particle of 
charge — e (e > 0) and mass m moves in the semi-infinite plane. A uniform magnetic field B 
is applied perpendicular to its surface. Cartesian coordinates are defined such that the x axis 
is perpendicular to the boundary and the motion is confined to positive values of x, y axis 
is along the boundary. Dirichlet boundary condition is imposed at x = 0. In order to work 
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with a finite system, we introduce another Dirichlet boundary condition at x = L±, where 



L± is taken so large with respect to the magnetic length Ib = yhc/eB that an electron in 
an eigenstate cannot feel the presence of both boundaries. In this case it is clear that the 
contributions of the two boundaries to g{s) will be identical (we will check this explicitly 
below). For the same reason of working with a finite system, we impose periodic boundary 
condition in the y direction, taking L to be the length of the boundary. 
In the Landau gauge A = (0 , Bx), the Hamiltonian of the particle is: 



The Green's functions G{E] r, r ' ) and Goo{E; r, r ' ) being periodic in y — y' one can expand 
them as Fourier series, e.g.: 

G{E- r, r ' ) = ^ G,^ {E-x,x'). (17) 

Py 

The sum runs over all Py = Uyh/L, where Uy G 7L. Introducing the dimensionless variables 
X = \/2x/Ib, xq = yflpy/rmolB, = muL'^ /2Trh and v = — ^ — |, (p^ transforms into 



^ - ^ - ^ (5 + 5o)'^ Gp^{E; X, x') = -^{A7r%)-h{x - x') 



(18) 



The solution of this equation such that G'pj^(£'; 0, x' ) = Gpy{E; L±,x' ) = can be written 
as a sum of three terms: 

2777, ~ ~ ~ 

Gp {E;x,x') = — {Gooii^,xo;x,x') + G'ib(z/, xq; x, x') + G2b(z/, xq; x, x')). (19) 

The Green's function for unbounded motion Goo{i^, xq; x, x') is given by: 

Goo{jy,xo;x,x) = - ^^J^ D^{xy - xq)D^{xq - x<), (20) 

V ZTT 



where x> = max{x, x'}, x< = min{x, x'} (we used the Wronskian W{ D^{u), Dy{—u) ) 



27r 



The two last terms in (|T9|) are respectively the contributions of the boundaries 
X = and x = L±_ to the Green's function. They are found by requiring that Gpy{E; x, x' ) 
have to vanish on the boundaries and neglecting exponentially small terms due to the 
asymptotic behaviour of the function as its argument is large and positive (recall that 
L± = V2L^/Ib > 1): 

Gib(z/, xo; X, x) = . D^[x - xq) 

V27r D^[-Xo) 

G2,{u, xo; i, i) = ^^^4^^^^D^(xo - X). (21) 
V27r L)^(xo - L±) 

Since we are interested in the limit L — » oo, we can now replace the sum over py in (|T^ by 
an integral. We deduce that Gib(i/, xq; x, x) and G2b{y, xq; x, x) give the same contribution to 
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(p!5|), which can be expressed using the relevant formulas |]T3| for the integrals of the function 

D„ as: 




-xq) - XoV^(- 



(22) 



This is the resolvent for the semi-infinite plane. In the infinitely large strip geometry, 
g{s) — Qoois) is twice this result. 

Formula (^) can be understood more intuitively as follows. The integrated density of 

states N{E) count the number of zeros of the wavefunction at the origin 1(2^0)) ior all 

Xq and all energies e up to e = -E: 

C{E) 

The sum and the integral cannot be inverted since there are infinitely many states. However, 
we can formally subtract from N{E) the integrated density of states of the infinite plane 
which we write as: 

N^(E) = \^^^(E-Mn.^\» = -^Z f ^'(-i;^\>- (24) 

Py Py ^^^^ 



equation (j^) we obtain: 



a factor | is introduced since we neglected the second boundary contribution in Using 



9is) - goois) = l-J2{-^\n D^{xo) - ) (25) 

Py 

Replacing here the sum over py by an integral we obtain Eq. (p^). The heat kernel is found 
by performing inverse Laplace transform on (|2^) or (|25|). 

It is interesting to note that g{^/E) and P{t) have the following forms: 

hujgiVE) = N^9oo{^) + ^/Nigsi^) 

P(t) = p^{ujt) + ^P5(cut), (26) 

where iV<i> = mujLL±/2TTh, Nl^ = muL'^ /2TTh and none of the functions decorated by tilde 
depend explicitly on B, L or Lj_. The boundary term is proportional to \jNj^ and smaller 
than the bulk term by a factor of order ;^^= for L± ~ L. For finite lengths L of the 



boundary, there are exponentially small correction terms in p^) and pq) due to the error 
introduced by replacing the sum in by an integral in (pT]). Neglecting these small corrections, 
an asymptotic expansion of the heat kernel and of the density of states in h^/f = l/iV$, 
which is a small parameter (assuming large fields or large systems), has two terms (this is 
due to the straight boundary - there will be more terms if an edge has a non-zero curvature). 

The arguments of this section extend to any separable system provided one can reduce 
the problem to a one-dimensional Sturm-Liouville problem. As an example we derive in the 
appendix the resolvents of a particle moving in a disc with and without a magnetic field. 
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IV. WEYL EXPANSIONS OF THE RESOLVENT AND THE HEAT KERNEL. 



In this section we shall derive asymptotic expressions for the heat kernel and the re- 
solvent using equation (p5|). For large positive energies (large negative z/) and large Xq 
the Darwin asymptotic expansion (basically WKB expansion of the wavefunction) of the 



parabolic cylinder function [14,13] is 



lnD,(xo) ~ ^ - ^ In T{-u) - 9{io, a) - ^H^l + 4a) + £ ^^-^^ , (27) 

where 9{xo,a) is an odd function of Xq, a = ^ > and Xq + 4a ^ 1. The coefficients 
are odd functions of Xq for odd values of s and even functions of Xq for even s. Since 
we integrate Xq in (^5]) over a symmetric interval, only even functions of Xq contribute. 
The first three even-indexed coefficients d^g are given respectively by [Q: d^ = |xq — 2a, 
d^2 = - 186a£2 + 80a2 and dis = ^x^ - 29862ax^ + 62292a2£2 _ 3i|32^3_ 
Using (j2^) and (pT]) the asymptotic expansion of the resolvent is obtained: 



/T — , N^, I 7 

q{s = vnuja) ~ - — H + 

' huo^ 24a2 960a4 



ttAT; 1 9 2625 241197 

+ a -i::;r^ + ■■■)■ (28) 



Ahu 256a5 262144a2 225a 

The Weyl expansion of the heat kernel is derived using (||) 

r iV^ r2 Tt^_ , ^5r^_ 7309r6 

^ J~ 24 + 5760 64+ 16384 315 220 + ' ' '^^ ^^^^ 

These expansions could in principle be continued indefinitely, calculating recursively the 
coefficients d^s- 

It is instructive to compare the expression ( ^9]) with its counterpart at i? = 0. The heat 
kernel in the absence of magnetic field is simply obtained from (|29| ) by taking the argument 
of the analytic function of r in each parenthesis to be zero. This suggests that a natural 
extension of the heat kernel expansion for the two-dimensional shape without the magnetic 
field given for example in [^,^ to the case of magnetic billiard is to multiply each term of 
that expansion by an analytic function of r. These functions should be universal for all flat 
two-dimensional billiards with smooth boundaries. 



V. THE BALIAN-BLOCH METHOD. 



The small time asymptotic expansion of the heat kernel can also be found using a method 
suggested by Balian and Bloch 0. It consists in a reformulation of the problem of solving a 
partial differential equation of elliptic type with Dirichlet (or another) boundary condition 
in terms of an integral equation of the Fredholm type. This integral equation is then solved 
iteratively (Neumann series), and each term in this multiple reflection expansion corresponds 
to one term in the asymptotic series of the heat kernel, as shown by Balian and Bloch |^. 
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The Green's function for the Dirichlet problem in a domain 5* with a boundary dS in 
absence of magnetic field is written as a sum of two terms: 



G{E; r,r') = G^E; r, r' ) + GsiE; r, r' ), (30) 



where each term in the right-hand side satisfy ([1^) in which H = — is the Laplacian. 
The first term is the infinite plane Green's function; the second is specified by the boundary 
condition: Gs{E]r,r' ) = —Goo{E;r,r') for r on dS. This boundary term is expressed in 
terms of an unknown density yU£;(Q;,r) as: 

Gs{E; r,r ) = J da^ — /i£;(Q:, r ), (31) 

as " 

and /i£;(Q;,r) is determined by solving the following Fredholm integral equation: 

^/i^(a,r') = -Goo(i?;a,r')-|rfa^^^2£^^^MAr'), (32) 

dS 

where q;,/3, . . . are arbitrary points on the boundary dS, dap is the boundary differential 
element, and is the normal derivative at the point (3 with the normal oriented towards 



the interior of the domain. Solving iteratively the integral equation (|32D for the density 
/i£;(Q;,r), the following multiple reflection expansion is obtained for the Green's function: 

G{E- r, r' ) = G^E; r, r' ) - - / daJ^^^^^^^^G^{E- a, r' ) + 

m J orir, 
as 

(g) ' I ,..,. /^-(^^ ■■■ ") ^) CUE; py)-... (33) 

QS np 

By mean of the time-dependent Green's function G^{t; r, r') introduced in section 4, we can 
calculate the heat kernel: P{t) = Jg d^vG^{t\ r, r). The multiple expansion of G'^{t\ r, r') is: 

G+(t; r, r' ) = (t; r, r' ) - ^ | | dr^-^^^^^p^G^r- a, r' ) + 



[m ) J ^^'^^^^ J J d^^ d^, ^' ^ ) - • • • (34) 

^ ^ as 

This approach was applied in this form to the case with a uniform magnetic field per- 
pendicular to the domain S j^. However, at each order of the obtained multiple reflection 
expansion we obtain a term which is not gauge invariant. This problem may be easily cor- 
rected by introducing the covariant derivative ^ ^An{ct) instead of the usual This 

substitution is of no importance when the gauge is chosen such that the vector potential has 
no component normal to the boundary (as it happens to be in our problem), but generally 
should be taken into account. 

Let us therefore apply the Balian and Bloch method to the semi-infinite plane in a 
uniform magnetic field. The Green's function for the infinite plane is given by: 
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(35) 



We calculate now the first term (proportional to B"^) in the small magnetic field expansion 
of the heat kernel. It turns out that the first three boundary-dependent terms in the multiple 
refiection expansion do contribute to this order. The small B expansion of the one-refiection 
term begins with: 

L / 2m L / 2m 7 22 

for the two-refiection and three-refiection terms, it is respectively: 

L / 2m 3 99 
8 V Trht 192 

and 

L 2m 1 99 

— \ — uH^. 

8 V ttM 192 

Therefore the small magnetic field expansion of the heat kernel begins as follows: 

L 2m L / 2m 899 , , 

P(t)=P^(t) \/^ + -\/^ cuV-... (36) 

in accordance with the previously obtained result (PD|). This however disagrees with Rob- 
nik's calculation 0. In this work the author assumed that in the case of a zero-curvature 
boundary, the one-refiection term contains all the correction due to the boundary (as it is for 
the problems with zero magnetic field). However as already noted by John and Suttorp [Q, 
higher order terms can not be neglected in calculating physical quantities in the geometry 
with straight boundaries in the presence of a magnetic field. We also remark that in order 
to calculate the term proportional to u'^"' in the asymptotic expansion of the heat kernel, all 
the multiple refiection terms up to the 2n + I'th order have to be taken into account. This 
result is intuitively appealing. Indeed as the magnetic field increases, the trajectories of the 
particles bend more and more, so that higher and higher terms in the multiple refiection 
expansion do contribute. 



VI. PERIMETER CORRECTIONS TO THE LANDAU DIAMAGNETISM. 

We apply in this section the results of the preceding sections to compute the magnetic 
susceptibility of an independent electron gas on a strip of very large width L±. 

Let us first consider the case of a non-degenerate gas at temperature T = The 
heat kernel P{t) for t = t/u = hjS is nothing but the one-electron partition function. If N 
is the electronic density (per unit area) and = 2^ Bohr magneton, the magnetic 
susceptibility of the ideal gas is given by 

X = \immfilP—\nZ. (37) 
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Using the expansion ( p9D of P{t) we obtain the weak field susceptibihty, 



where Xt = y ^^^^ is the de Broghe thermal length. Like in the correction to the Landau 
diamagnetic susceptibility Xoo = ~\Nl3ii\ is paramagnetic, but is smaller by one order of 



magnitude. Our result agrees with numerical calculations [|T^. 

For a degenerate electron gas of Fermi energy Ep-, the connection between the suscepti- 
bility and the partition function is different 0. In the grand-canonical ensemble, x is given 
by 

VL{Ef) being the thermodynamic potential at the chemical potential Bp- The thermody- 
namic potential at zero temperature VLq{Bf) is connected to the partition function (heat 
kernel) through p 

joo+0+ 

^.{E,) = - I ^^e^^-p-'ZiP). (40) 

-ioo +0+ 



Using ( [26| ) we obtain: 



- -If T l^i' - ("^-^^-<^) - S--'^^'^)) ■ 

-ioo+0+ ^ V -L / 

Replacing in this equation the Weyl expansions (|29|) of Zooij) and 2'5(r) amounts to ne- 
glecting all oscillating contributions due to periodic orbits in VLq{Bp). However, even if 
the oscillating terms in the thermodynamic potential can at very low temperature give the 
largest corrections to the Landau diamagnetism, these terms are exponentially damped by 
temperature IJl^,!!!] (more precisely, the contribution of each periodic orbit is damped by 



the factor ^p^^ sinh ^(^p^^); where r is the number of repetitions of the orbit and 

kp = ^^^"^^ is the Fermi wavevector). For temperatures such that -C ksT <^ Bp, one 
can neglect these oscillating contributions, and we have at small field: 

= ~247rmc2 " wTTk^) ' ^^^^ 



Here again the correction is paramagnetic and coincides with perturbative calculations 

The expressions (^) and (|4^) give the magnetic susceptibilities of a semi-infinite plane 
regularized to have a finite width L_i_ for small fields B such that L_i_ ^ Ib- Here again the 
second term in these equations have to be multiplied by two if one wishes to consider a finite 
very large strip (with two boundaries) instead of the semi-infinite plane. 

The susceptibilities (|38|) and ( ^2]) give the perimeter corrections to the Landau diamag- 
netism for a general macroscopic billiard of surface 5* and smooth boundary of length L, 
with Lj_ = S/L. However for a generic billiard this correction is not the only one, and other 
terms due to the curvature of the boundary must be added. 
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VII. CONCLUSION 



We obtained an exact relation for the resolvent and the heat kernel in the semi-infinite 
plane geometry in the uniform magnetic field. This relation is a convenient starting point 
for the asymptotic expansion, using the WKB approximation for the wavef unctions. The 
validity of the method extends beyond the special case of the semi-infinite geometry, in 
fact it works for any separable system. The main result of the work is the boundary term 
of the heat kernel asymptotic expansion, which is an infinite series calculable recursively 
in principle. The first terms in this series were verified using the Balian-Bloch multiple 
scattering expansion, which is less convenient practically than the method of Stewartson 
and Waechter. We point out how the presence of the magnetic field affects the physics of 
the Balian and Bloch approach making necessary to consider three terms in the multiple 
scattering series in order to get correctly the weak field response. The corrected perimeter 
corrections to the Landau diamagnetism are derived. 

The properties of the heat kernel of manifolds have been extensively investigated starting 
(among others) from the work of Kac in order to relate the spectral and geometrical 
descriptions [^. Our results could be seen as an extension of these works to the case of 
magnetic billiards. For the case of a straight boundary, the magnetic field adds an infinite 
series to the bare perimeter term, which is, in some sense, equivalent to an effective curvature 
of the boundary. If the boundary has a curvature, another length scale enters the problem, 
coupling with the cyclotron radius. 
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A. APPENDIX 



Consider a particle confined to a disc of radius R. Introducing the momentum k"^ = ^p^, 
the energies are solutions of the equation Ji{kR) = 0, where I is the angular momentum 
quantum number. The number of zeros is equal to the integral of the logarithmic derivative 
of the Bessel function: 

N{E) = Y. f ^±4M. (A.l) 

' C{E) ^ ' 

Then according to our method, we obtain the resolvent 

^^^^ -\irk -um- ^ ^ ^ 

which gives the result of Stewartson and Waechter |^ after subtracting the resolvent corre- 
sponding to the infinite plane. 

In the presence of a magnetic field, the energy of the particle moving in a circular billiard 
is given by: — 1 + |/|; A$) = 0, where A$ = The counting function is now: 
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' C{E) 2 nw' ^ I'h 



and the resolvent is: 

We emphasize again that in order to work with well-defined quantities we must subtract 
from this expression the part of the resolvent corresponding to the infinite plane. 
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